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Abstract. We show that the extinction of quasar absorbers increases exponentially with the 
logarithmic column density of any volatile metal (e.g. zinc), with a charateristic turning point 
above which the quasars are severely dimmed. We derive the relation between extinction, HI 
column density (A''hi), metallicity {Z = Nznii/Nm) and fraction of iron in dust [/fo(Z)] in 
Damped Lyman a (DLA) systems. We use this relation to estimate the effect of dust obscuration 
on the statistical distributions of A'hi and Z measured in magnitude-limited surveys of DLAs. 
In the redshift range where the measurements of zinc column densities have sufficient statistics 
(1.8 ^ z 3) we find that the obscuration bias affects the shapes of the observed distributions. 
The metallicity distribution is particularly affected by the bias, which hampers the detection of 
DLA galaxies with near solar metallicity. Our results perfectly reproduce, without tuning the 
dust parameters, the DLA observational threshold logA'^(Znll) [cm~^] < 13.15 found by Boisse 
and collaborators in 1998, in terms of a rapid rise of the obscuration. Our predictions of the 
effects of the bias on the statistics of DLAs are consistent with observational results obtained 
from unbiased surveys of radio-selected quasars. 
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1. Introduction 

A variety of astronomical observations indicate that interstellar dust is a pervasive 
galactic component, not only in the Milky Way and nearby galaxies, but also in the high- 
redshift Universe (Meurer 2004). Among quasar absorbers. Damped Lyman a (DLA) 
systems are the best candidates to host dust because of their high HI column density, 
which shields the grain from ionizing photons, and of their metal content, higher than 
in other absorbers. If dust is present in DLA systems, it will scatter and absorb the 
radiation of the background quasar, dimming its apparent magnitude (extinction) and 
changing the slope of its energy distribution (reddening). In the most extreme cases, the 
extinction may obscure the quasar leading to a selection effect, the obscuration bias, 
which may prevent the detection of metal-rich, dusty DLAs (Ostriker & Heisler 1984, 
Fall & Pei 1993). 

Measuring the extinction or reddening of DLAs is a difficult task owing to the uncer- 
tain, variable continuum of the quasars. Detection of quasar reddening from intervening 
DLAs was first claimed by Pei et al. (1991). Recent investigations based on large data 
sets find reddening from metal absorbers at redshift 1 ^ z ^ 1.9 (Khare et al. 2005), but 
not from DLAs at 2 ^ 2.8 (Murphy & Liske 2004). Differential reddening in two lines of 
sight of one DLA a.t z — 0.9 has been reported from the study of a gravitationally lensed 
quasar (Lopez 2005). 

Evidence of the obscuration bias has been searched by comparing the statistics of 
magnitude- limited surveys with those of complete surveys of DLAs lying in front of radio- 
selected quasars. The CORALS survey has provided a marginal signal of obscuration 
in the absorber number density, the neutral gas content (Ellison et al. 2001) and the 
metallicity (Ellison 2005), but only at « la level in each case. 
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The lack of DLA systems with high values of metallicity and HI column density, above 
the threshold logA^(Znll) [cm~^] ~ 13.15, was tentatively attributed to the effect of 
dust obscuration (Boisse et al. 1998). The original motivation of the present work was 
to understand if that threshold can be quantitatively explained in terms of a fast rise 
of the DLA extinction, in which case it would represent an important clue of the obscu- 
ration effect. With this goal in mind, we first derived the relation between interstellar 
extinction and A^(Znll). We then used this relation to quantify the effect of the bias on 
the distributions of column densities and metallicities of DLAs. Finally, we applied our 
method to the current sample DLAs with Znll measurements. These three parts of our 
work are briefly outlined in the next three sections. More details can be found in the full 
paper (Vladilo & Peroux 2005). 



2. The extinction law 

The extinction in a photometric band with effective wavelength A, A\, is related to 
the line-of-sight zinc column density, Nzn = iV(Znll), by means of the relation 

= aAlO'°sJVz„ (2.1) 



where 

aA = 1.817 X 10-24 Ape X ffA X . . x(^\ . (2.2) 

(l-/zn) VZny 

These relations can be derived by comparing the column densities of a refractory and 
a volatile element in a given line of sight (we adopted iron and zinc, respectively). The 
expression (|2.2|l is composed of: (1) constant factors (e.g. Apo is the atomic mass of 
iron); (2) a dust grains parameter, gx; (3) the fraction of iron and zinc atoms which are 
incorporated in the dust, /pe and fzn', (4) the relative abundance of iron and zinc in the 
medium, (Fe/Zn). This last term is expected to be constant, close to the solar value, in 
the local ISM. The term /fo/(1 — /zn) is approximately constant and « 1 by definition 
of refractory and volatile elements. The only factor expected to vary is 

5A = s , (2.3) 

Tgr Pgr 

where s is a geometrical factor (=3/4 for spherical grains), Qe,\ the extinction efficiency 
factor (Spitzer 1978), rgr the grain radius, Qgr the grain density, and Xp^ the abundance 
by mass of iron in the dust. 

We tested and calibrated the relation H2.1|l using interstellar measurements of extinc- 
tion and A^(Znll). The results are shown in Fig.Q] The data follow the expected trend 
and yield a mean value of the ax parameter in V band ay — 0.29(±0. 07)10"^^ magni- 
tudes cm~2. From the condition dA\/d{logNzn) — 1 magnitude, we obtain a turning 
point log A^(Znll) ~ 13.18, above which the extinction experiences a dramatic rise. This 
extinction barrier lies remarkably close to the Boisse threshold and this fact motivated 
us to study the extinction relation for DLAs. As a result, we obtained that the observer's 
frame extinction of a DLA at redshift z is 

Ax{Nh, Z; z)^AoG e(^YT^) ho{Z) iVn Z , (2.4) 

where Ao ~ 1.85 x 10"" mag cm^, G = (.gA)DLA/(5A)iSM, and ^(A) = A{X)/A{Xv) 
is a normalized extinction curve; Z = Zn/H and A^h are the metallicity and HI column 
density of the absorber; fpdZ) is an analytical function, based on depletion data (Vladilo 
2004), that gives the fraction of iron in dust as a function of the metallicity. Since fpciZ) 
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Figure 1. Left panel: V band extinction versus zinc column density in the interstellar medium 
(see Vladilo & Peroux 2005). Right panel: Fraction of iron in dust versus metallicity in DLAs 
(Vladilo 2004); smooth line: analytical approximation adopted in the present study. 



vanishes at low metallicities (Fig. 1, right panel), the predicted extinction will accordingly 
vanish. Thanks to this fact, we do not need to know the value of the grain parameter 
at the very early stages of chemical evolution, when the behaviour of the dust is most 
unpredictable. In the application of Eq. 1)2. 4|l we consider both a 'Milky- Way' and a Small 
Magellanic Cloud ('SMC') type of dust. For each type we adopt an appropriate ^(A) and 
G, making also use of SMC extinction data (see details in the full paper). 



3. Method 

At a given redshift, the distribution of DLA extinctions and the obscuration bias are 
completely determined from l|2.4|l for any assigned pair of iV(HI) and Z distributions. We 
used this property to find the relation between true and biased distributions of iV(HI) 
and Z in DLAs. Simple mathematical expressions can be derived by assuming that the 
distributions of iV(HI) and Z are statistically independent and by ignoring multiple DLA 
absorptions in a given line of sight. Both assumptions are conservative in the sense that 
they may slightly underestimate the effect of quasar obscuration. The fraction of DLA- 
quasar pairs that can be detected in the differential element dNu dZ dz is 

n(m; z) dm 

B^,iN^,Z,z) = ^ , (3.1) 

n{m; z) am 

where A\{N}i, Z, z) is the extinction (|2.4(l , n{m) is the distribution of the quasar apparent 
magnitudes m, and mi the limiting magnitude of the survey. In the redshift interval 
zi ^ z ^ 22 the relations between the true distributions of iV(III) and Z in DLAs, /atjj 
and fz, and the corresponding biased distributios, /^^ and are 

^ J^-B„,,{N^,Z)hdZ 

Jo 



,b lZ^,B^,{Nn,Z)fN^dNH 



Jz — Tod tt; JZ , (.J-OJ 

where BmiiNn, Z) ~ Bm^iNii, Z,z) and TVdla = 10^°-^ atoms cm^^. The goal of our 
procedure is to invert these equations in order to recover the true distributions given 
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Figure 2. Flow diagram of tiie procedure adopted witii tlie aim of recovering tlie true distribu- 
tions of HI column densities and metallicities in DLAs starting from the observed distributions. 



the biased ones, which can be estimated from niagnitude-hmited surveys. To do so we 
use the procedure sketched in Fig. 2. We model the unknown distributions /jvh ^J^d fz 
using analytical functions with free parameters. We start from a set of trial parame- 
ters and compute the extinction distribution, the function B„h{Nyi, Z) and the biased 
distributions. We then compute the deviation between the biased distributions and 
the observed ones. The procedure is iterated to find a minimum in the values. At 
each iteration we estimate the true distribution of quasar apparent magnitudes, n(m), 
by correcting the observed magnitude distribution taking into account the effect of the 
extinction. In this way we do not need to assume a quasar luminosity function. For this 
and other reasons our procedure is different from the one presented by Fall & Pei (1993). 



4. Results 

As a first application of our method we considered the sample of DLAs with Znll data. 
The distribution of apparent magnitudes was measured from the quasar counts in the 
r and g bands of the Sloan Digital Sky Survey. The extinctions were estimated in the 
same two bands. The selected sample is characterized by a quasar Hmiting magnitude 
m£ ~ 19 in the optical, required to detect the Znll Hues at high spectra resolution. 
Most data are concentrated in the range 1.8 ^ z ^ 3.0, with a median redshift z « 
2.3. The current limitation of the statistics (about 40 DLAs in 2004) allowed us to 
compute the empirical distributions in a small number of bins. Because of the limited 
empirical contraints we modeled the true distributions Jn-h and fz using functions with 
the smallest possible number of parameters. For the A^(HI) distribution we adopted 
a power law /at^ = C x N^'^ truncated at the highest observed values of DLA column 
density, log (HI) ~ 22. For the metallicity distribution we adopted a Schechter function 
fz~Cx {Z/Z„)"e^^/^' . With tighter observational constraints available in the future, 
it will be possible to model the true distributions more realistically. 

For the different types of dust (MW- and SMC-type) and photometric bands consid- 
ered, we found that the slope of the A^(HI) power law lies in the range 1-5 < /? < 1.6, with 
typical fit errors of ±0.12. The mean value < Z > of the Schechter metallicity distribution 
lies in the range —0.4 < log(< Z > /Zq) < —0.3, with fit errors of ±0.15 dex and total 
errors (including the uncertainty of the empirical distributions) of ±0.3 dex. 



Obscuration of DLAs 



5 




0.0 



[Zn/H] 



Power law (^=1.61) 
M'-..... i.r.r; Schechter <a=-0.46 knee=-0.19) 

.iUi i.>p -Msl; r band «xUnction 



30 



3,0 0,2 0.4 0.6 0.8 1.0 

Extinction (magnitudes) 



Figure 3. Left: distribution of metallicities in DLAs for different predictions of our models (see 
Vladilo & Peroux 2005); the error bars indicate, from left to right, the ilcr errors of the mean 
metallicity of DLAs estimated (1) at face value, (2) from the CORALS survey (Ellison 2005), 
and (3) from our predictions. Right: distribution of DLA extinctions in the r band obtained for 
the Znll sample of DLAs for an SMC- type dust; vertical bar: median extinction. 



The mean metallicity that we derive is higher than the canonical value of DLAs taken 
at face value, < [Zn/H] >~ —1.1 dex. As we show in Fig. 3, our estimate is consistent at 
PS la level with the mean metallicity of the unbiased CORALS sample (Ellison 2005). 
In the full paper we show that also our predictions of the number density, ?t.dla(-z), and 
gas content, f^DLA, are consistent with CORALS results (Ellison et al. 2001). 

An example of distribution of DLA extinctions predicted by our study is shown in Fig. 
3 (right panel). The distribution has a median value of 0.14 mag in the r band, but is 
extremely asymmetric. Most DLAs have negligible or little extinction, but the long tail of 
the distribution, extending up to very high values, indicates that a significant fraction of 
DLAs may indeed obscure the background quasar. We estimate that the missed fraction 
is ~ 30% to 50% at limiting magnitude = 19, but decreases steadily with m^. Surveys 
mi > 21 are weakly affected by the bias. 

The extinction distribution highlights the difficulty of measuring the 'mean' DLA red- 
dening using large data sets of quasars. The predicted median reddening is relatively low 
(~ 0.05 magnitudes) and hard to detect on the uncertain quasar continuum. In addition, 
the procedure of stacking a large number of spectra for detecting the 'mean' redden- 
ing can give a meaningful answer only if the mean change of continuum slope is a fair 
representation of the change of slope of individual cases. 

The metallicity dependence of the iron dust fraction (Fig. 1, right panel) indicates 
that DLA extinction should become negligible at [Zn/H] < —1.5. In conjunction with 
the trend of decreasing [Zn/H] with redshift (Vladilo et al. 2000), this in turn suggests 
that DLA extinction may vanish at z > 3. Therefore, the difficulty of detecting quasar 
reddening from DLAs at z « 2.8 (Murphy & Liske 2004) may be due, in part, to the 
intrinsic decline of the dust content at high redshift. 

An interesting prediction of our method is the fraction of obscured DLAs as a function 
of A^(Znll), (j){Nzn)- In Fig. 01 (left panel) one can see that (t){Nzn) shows a sharp rise ex- 
actly around the Boisse's threshold log iV(Znll) ~ 13.15. This remarkable result, obtained 
without tuning the dust parameters, demonstrates that the obscuration effect provides 
a natural explanation for the threshold. To discard this explanation one should invoke 
some unknown mechanism for explaning the lack of DLAs with log A^(Znll) > 13.15 and, 
in addition, one should tune the dust parameters in such a way to move the threshold at 
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Figure 4. Predicted fraction of obscured DLAs as a function of total Zn column density (left 
panel) and metallicity (right panel) in a DLA-quasar survey with limiting magnitude ~ 19 for 
different models (MW- and SMC-type dust; r and g photometric bands). 



much higher values of log A'^(Znll). Our explanation is certainly more attractive in terms 
of economy of hypothesis. 

The reality of Boisse's threshold has been invoked to bring into agreement predictions 
of galactic chemical evolution models (Prantzos & Boissier 2000) and of cosmological sim- 
ulations (Cen et al. 2003) with seemingly discrepant observations of DLAs. The existence 
of a physical explanation for the threshold favours a scenario in which the predictions of 
such models are generally correct and the obscuration bias plays an important role. 

The combination of the hmit logA^(Znll) < 13.15 with the cutoff logA^(HI) ^ 20.3 
imposed in the surveys creates a cutoff of the highest possible metallicity detectable in 
DLAs. The way the obscuration bias affects the metallicity distribution is shown in the 
right panel of Fig. 0] One can see that the bias severely hampers the detection of DLAs 
with near solar metallicity. Because of the implications on our understanding of the star 
formation history of the Universe, the effect on the metallicity distribution is probably 
the most important consequence of the obscuration bias. 
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